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m ■ Abstract 
> 

t^J- ■ Modifications of the fi + lifetime in matter due to muonium (M = /i + e ) 

formation and other medium effects are examined. Muonium and free u + 
qs, , decay spectra are found to differ at 0(am e /m^) from Doppler broadening 

and 0(a 2 m e /mn) from the Coulomb bound state potential. However, both 
& types of corrections are shown to cancel in the total decay rate due to Lorentz 

and gauge invariance respectively, leaving a very small time dilation lifetime 
^h. difference, (tm — t ^ l +)/'t^+ = °2mZ e — 6 x 10~ 10 , as the dominant bound state 

effect. It is argued that other medium effects on the stopped /u + lifetime are 
similarly suppressed. 

The muon lifetime, r„, is very well measured. Its current world average 



= 2.197035(40) x 10" 6 sec (1) 

exhibits an uncertainty of only 18 ppm. From that lifetime, the Fermi constant (denoted 
here by G^) is determined via the defining relationship 

f(x) = 1 -8x + 8x 3 -x 4 - 12x 2 lnx, (2) 

where m M is the muon pole mass and R.C. stands for QED radiative corrections to muon 
decay as calculated in an effective local V-A theory. Other standard model electroweak loop 
corrections as well as possible "new physics" effects are absorbed in G^. 

The R.C. in (||) have been computed [fJH through 0(a 2 ) and higher order logs have 
been obtained using the renormalization group [|J. Altogether, one finds [^,|J 
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7r V3 m e J ' 7r 2 V9 m e m e 
a' 1 = 137.03599959(40) (3) 

where C corresponds to unknown non-logarithmic 0(a 3 /7T 3 ) ~ 10~ 8 corrections which are 
assumed to be insignificant. Employing ([|,@,^D leads to 

= 1.16637(1) x 10~ 5 GeV" 2 . (4) 

The Fermi constant can be compared with other precise measurements such as a, mz, 
sm 2 9w, ttiw, etc., and used to test the consistency of the Standard Model at the quantum 
loop level. For example, G M = ixaj \[2m\ v {\ — mfp / 'm|)(l — Ar) where Ar ~ 0.0358 (for 
m top — 174.3 GeV and Higgs mass niu — 125 GeV) represents calculable electroweak radia- 
tive corrections 0. In that role, G M helped predict the top quark's mass before its discovery 
and currently constrains the Higgs mass to relatively low preferred values m# < 220 GeV. 
In addition, it provides a sensitive probe of "new physics" such as SUSY, Technicolor, Extra 
Dimensions etc. 

Recently, there have been several proposals |8|-[l0l to further improve the measurement 



of (and thereby G^) by as much as a factor of 20, bringing its uncertainty down to an 
incredible ±1 ppm. Of course, such an improvement can only by fully utilized if the other 
electroweak parameters with which it is compared reach a similar level of precision and 
radiative corrections to their relationship (e.g. Ar) are computed at least through 2 loops. 
That confluence of advances appears unlikely in the foreseeable future. Nevertheless, given 
the fundamental nature and importance of G^, efforts to improve its determination should 
be strongly encouraged and pushed as far as possible. In that spirit, we examine in this 
paper several theoretical concerns that must be addressed in any ±1 ppm study of r M . 

Precision muon lifetime experiments involve stopping /i + in material. At some level, the 
medium in which it comes to rest will affect the muon decay rate. The most straightforward 
issue to consider is the formation of muonium (M = /i + e - ) and its bound state effect on r M . 

During the slowing down of fi + in matter, muonium will form and be ionized many times 



il"|,[T2"|| . Eventually it comes to (thermal) rest at which point the muon or muonium atom 



has a room temperature kinetic energy ~ 0.04 eV or thermal velocity /3* hermal ~ 2.7 x 10~ 5 



The actual fraction of stopped /i + that end up as muonium and decay while in that 
bound state is very medium dependent. It can range from a small fraction in metals to 
nearly 100% in some materials. If a significant modification of r M occurred in muonium, 



a correction would have to be applied. For example, a published study |T3[ has claimed a 
0.999516 reduction factor for the bound state decay rate. Such a large —484 ppm shift would 
be difficult to correct for at the ±1 ppm level unless the muonium formation fraction was 
very precisely known. It would also impact the interpretation of existing measurements 
fl~4]| and might indicate the possibility of other large medium effects. 

Given its potential importance, we have reexamined the muonium bound state effect on 
the muon lifetime. As we shall show, the leading correction turns out to be remarkably 
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small, ^y?r — 6 x 10~ 10 and easy to calculate. 



We begin by recalling some basic properties of muonium [|T^]. It is a fi + e Coulombic 
bound state with hydrogenic features. The reduced mass 

m e m,, , , 

m = ~ 0.995m e 5 

m M + m e 

is slightly below the electron mass, but the difference is inconsequential for the considerations 
here and can be neglected. The energy levels are given by 

a 2 m , n 

£n = -y^- = -13.5eV/n 2 , n=l,2... (6) 

For the n = 1 ground state the virial theorem tells us that the average (electron) kinetic 
energy is 

(T) — —Ex — 13.5 eV (7) 
while the average bound state potential energy is 

(V) =E X - (T) = a 2 m = -27.0 eV. (8) 
Also, the n — 1 momentum distribution is given by |TJ| 



Aitp 2 \^(p)\ 2 = —am{ - am 9 - ) p 2 . (9) 



7r \p 2 + a z m 

It peaks at p 2 = a 2 m 2 /3 which is somewhat below the average (p 2 ) = a 2 m 2 . In that 
configuration the e" and fi + have equal but opposite momenta and r.m.s. velocities, /3 r . m . s . = 
(/3 2 ) 1 / 2 , f3 e = a ~ 1/137, and (3^ ~ am e /m^ ~ 3.5 x 10~ 5 . Note that the muon's r.m.s. 
bound state velocity is similar (somewhat larger) in magnitude to its room temperature 
thermal velocity. 

Of course, (^) represents a non-relativistic approximation and should not be used indis- 
criminately for large p. However, it is adequate for the analysis presented here. In fig. [I], we 
display the bound state momentum probability distribution corresponding to (0). 

Muonium bound state modifications of the muon decay spectrum and lifetime must 
vanish as a or m e /m M — > 0. It is, therefore, useful to organize such corrections as an 
a n (m e /m^) m , n,m = 1,2... expansion in the small parameters a ~ 1/137 and m e /m^ ~ 
1/207. Terms of order n + m > 4 are completely negligible and can be safely ignored. 

The leading 0{am e /m IJ ) effects result from the muon's non-zero velocity in the muonium 
rest frame. (For this discussion we take the muonium atom to be at rest in the lab frame.) 
The ground state velocity, (3^ ~ am e /m M , will modify the e + decay spectrum relative to its 
shape for free /i + decay at rest (see fig. |]). For a given it gives rise to Doppler broadening 
which corresponds to an o verall d ilation along with smearing of the positron energy over 
the range (E e + ± /3 fl p e +)/ — f3 2 . That distortion amounts to about ±1.9 keV near the 
free decay endpoint E e + = (m 2 + m 2 )/2m At for (3^ ~ am e /m^. Such an effect leads to a 
small spectral tail beyond the usual endpoint. Of course, that enhancement is accompanied 
by a depletion below the usual endpoint. In fact, smearing corresponds primarily to a 
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redistribution of the spectrum. Lorentz invariance requires that such effects cancel in the 
total integrated decay rate up to an overall time dilation factor 7 = 1/^1 — which slightly 
reduces the decay rate by I/7. Averaging over the momentum distribution in Eq. (§) or its 
relativistic generalization replaces (3 2 by (f3 2 ) ~ a 2 m 2 /m 2 and one finds |17| 



T(M -> e~e + u e ^) = (l- T(p + -> e+z/ e ^), (10) 

Although the bound state velocity effects on the spectrum are relatively large 0[am e jm p ) ~ 



35ppm compared to the ±1 ppm experimental lifetime goal, the overall shift in Eq. ( 12|) is 
completely negligible. So, lifetime counting experiments that are independent of spectrum 
shape details are very insensitive to small muon bound state velocity effects. 

In the case of muonium stopped in matter, the bound state and thermal muonium 
velocities should be approximately added in quadrature 

total \ 2 ^ ( oM\ 2 1 / /othermalA 2 



(0?) +(/?r mai ) • (13) 



At room temperature, / z3* hermal ~ 0.8/3^ and the total time dilation shift in Eq. (12) increases 



by about a factor of 1.6. Nevertheless, it remains negligible at the 1 ppm level. 

The muonium bound state Coulomb potential, (V) ~ — a 2 m e , also modifies the positron 
decay spectrum. Those 0(a 2 m e /m^) corrections correspond to ~ —27 eV modifications of 
the spectrum. If they were to also affect the muon lifetime at that order, the effect would be 
larger than time dilation and possibly near the ±1 ppm goal of future experiments. However, 
as we show in the Appendices, the leading effects cancel, leaving behind contributions that 
are much smaller than 1 ppm. 

Aside from Doppler broadening, the bound state Coulomb interaction has two main 
effects on muon decay: phase-space suppression, because the muon is off mass shell, and 
final-state e + e~ interactions after the decay. In fact these two effects largely cancel, because 
of charge conservation. We present a rigorous proof of this cancelation in the Appendices, 
but the nature of the cancelation is easily understood if we replace muonium by a simpler 
system, where the muon is bound to a static charge distribution with density p(r) = e\ip{r) | 2 , 
rather than to an electron. 

Initially the muon in our simplified system is at the center of the charge distribution 
with negligible kinetic energy and potential energy — Vq, where Vq = 27 eV. Since the muon 
is below mass shell, the maximum energy allowed to the positron after the muon decays is 
lowered by 27 eV from what it is for free-muon decay. This phase-space contraction tends 
to lower the decay rate. The positron, however, is formed at the center of our "muonium" 
atom and, because it has the same charge as the muon, it initially has the same potential 
energy, — Vb, by virtue of its interactions with the charge distribution p(r). The positron 
rapidly leaves the atom, but its asymptotic kinetic energy is reduced from its initial kinetic 
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energy by the 27 eV needed for it to climb out of the potential well due to p(r). Thus the 
entire positron energy distribution is shifted down by V = 27 eV: 



ar M (E) dr,(E + v ) 

dE e+ dE e+ { ' 

This is the effect of final-state interactions. The two effects, phase-space contraction and 
final-state interactions, cancel in the total decay rate because the energy distribution and 
the kinematic limit are shifted by the same amount, leading to the same total integral for 
the decay rate: 



W2-V0 dY M (E) dE = p/2-vb dT ^E + V Q ) ^ 



dE e + J dE P + 



At very low-energies, the positron can bind to the charge distribution and our analysis must 
be modified. However any such modification is of order o?m\jvr? <C a 2 m1/m'j l (the time 
dilation effect) for the total decay rate since the spectrum is highly suppressed, by phase 
space ~ pEdE, in that low-energy region (see fig. 2). Such corrections are negligible. 

Other effects, e.g. hyperfine interactions, e + e~ hard scattering etc., are also of higher 
order in a and/or m e /m^ and therefore suppressed. Thus, overall, the time dilation effects 
in (|10D— fljjD represent the dominant bound state correction to the muonium lifetime. 

The cancelation of 0(a 2 m e /m IJL ) effects in the total /i + decay rate is analogous to other 
cancelations of a similar nature that have been previously pointed out. 

A recent example of that phenomenon has been encountered in inclusive B decays. The 
total decay rate of a meson or baryon containing a heavy b quark can be approximated by a 
"free" b quark decay rate. If the rate is expressed in terms of a short distance b quark mass, 
corrections to that approximation are 0(1/ ml) rather than 0(1/ mi). They stem from time 
dilation, just as we have found for the much simpler muonium bound state, as well as from 
hyperfine interactions (which for muonium are suppressed by a 4 ). In that case, the absence 
of l/rtib corrections is due to color conservation ||18|| . Our NRQED analysis in the Appendix 
is closely related to the standard B analysis. 

Another example of the cancelation of V dependence in total decay rates was found by 
H. Uberall |L9|| in his classic study of muonic atoms, jJT— Nuclei bound states. There, the 



/i~ is bound to the nucleus by a much stronger V ~ —Z a m^ than muonium. Nevertheless, 
Uberall found that the Coulomb potential reduction of initial state phase space was canceled 
by final state e~— Nuclei Coulombic interactions in the total decay rate. In fact, he also found 
for that example that the leading bound state effect was time dilation which reduces the 
effective \i~ decay rate in muonic atoms by 

7" 1 = Vl - Z 2 a 2 ~ 1 - \z 2 a 2 } (Z 2 a 2 < 1). (16) 

Of course, besides ordinary decay \x~ — > e~u e u^, the muonic atom can undergo weak capture, 
H~p — * u^n, with protons in the nucleus. In fact, for heavy nuclei with Z > 12 the capture 
process dominates over ordinary decay. Capture rates are generally obtained by measuring 
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the fi lifetime in a material and subtracting out the ordinary decay rate (/i — > e z/ e *7t) 
using the correction in (16). That prescription needs some adjustment for high Z nuclei 



where residual low energy e —Nuclei interactions can significantly suppress the decay rate 



and become comparable to time dilation |20 



The analog of fi p capture for muonium is annihilation M — > VeV^. We have computed 



that rate for the n = 1 ground state and find [f2l] ,p2 



/ \ 3 

T(M -> v e v^) = 48vr e - \ T(/i + -> e + ^ e ^) 

~ 6.6 x 1(T 12 T(/i + -> all). (17) 
Our result is in accord with the more general analysis of Ref. That rate is of higher 



order in am e /m M than the time dilation correction in (|TT|) and 100 times smaller. It can be 
safely neglected. 

The very small [i~p — > v^n capture rate for protons has been determined by compar- 



ing measurements of r M + and r M - in liquid hydrogen [ 24"| . In that way, the induced weak 
pseudoscalar coupling g p of the nucleon can be inferred. Our general analysis along with 
Uberall's study lend theoretical support to this method. Increasing the precision of those 
measurements by a factor of 10-20 (if possible) could provide an interesting confrontation 
with theory |^5[ . 

Our analysis of the muonium bound state effect on the fi + lifetime provides a simple 
way of estimating other effects of a stopping medium. For example, in metals the fraction 
of stopped fi + in muonium is tiny. Instead, the thermal muon is screened by conduction 
band electrons. The density of screening electrons near the /i + is expected to be comparable 
to muonium, resulting in a similar V and average fi + velocity. Hence, the distortion of the 
e + emission spectrum due to Doppler smearing and the final e + potential energy should be 
roughly the same as muonium. Due to Lorentz and gauge invariance, the leading effects will 
continue to cancel in the total decay rate. The main medium effect on the /i + lifetime will 
again be time dilation due to (3^ ^ (from thermal and electromagnetic interactions) which 
will be negligible ~ 0(1O~ 9 ). In the Appendix, we describe how such medium effects can 
be analyzed using nonrelativistic QED. 

A possibility that we need only briefly discuss is the effect of the medium on radiative 
muon decay. Ordinary free muon decay bremsstrahlung, /i + — > e + v e v^ (as well as 77 and 
e + e~ production), is included in the 0(a) and 0(a 2 ) terms of Eq. (^). Those effects are 
rather small, despite the fact that soft and collinear bremsstrahlung can significantly modify 
the spectrum. As in our above analysis, those large logarithmic infrared effects primarily 
give rise to a redistribution of the spectrum but tend to cancel in the total decay rate p6| , |27| . 

Our general arguments carry over to radiative muonium decay. The radiative muonium 
decay spectrum will be distorted by Doppler smearing and V dependence, but those effects 
will largely cancel in the total decay rate, again leaving time dilation as the primary cor- 
rection. A similar argument can be applied to alternative medium potential effects such as 
screening in metals. 

Other long distance properties of a stopping medium may be more subtle to analyze. 
For example, after the /i + decay process, an outgoing e + might emit outer bremsstrahlung 
or Cerenkov radiation due to its material environment. Those effects will modify the e + 



6 



and electromagnetic emission spectra. However, such changes again correspond to a redis- 
tribution of events and should cancel (to a very good approximation) in the total decay 
rate or lifetime (assuming Lorentz and gauge invariance are maintained in the analysis). 
The space-time scale associated with the actual decay process (not the lifetime) is 
That space-time interval is much too small to be affected by long distance environmental 
conditions. 

In summary we have found that lifetimes of muonium and a free n + differ in leading order 
only by time dilation effects which are ~ 10~ 9 and completely negligible for future ±1 ppm 
experiments. The smallness of that difference follows from Lorentz and gauge invariance 
which guarantee cancelations among considerably larger spectrum distortions. A similar 
suppression of other medium effects on the fi + lifetime is also expected. 
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APPENDIX A: BETHE-SALPETER AND NRQED ANALYZES 



In the text, we described how spectral shifts occur in muonium decay but largely cancel 
in the total decay rate, leaving time dilation as the primary bound state correction. Here, a 
proof of those cancelations is given using a Bethe-Salpeter bound state approach. A second 
formalism, nonrelativistic QED (NRQED) [EH] , is then employed to illustrate the important 
role of electromagnetic gauge invariance or charge conservation in the cancelation. That 
method provides a powerful means of parameterizing other more general medium effects 
and showing that they also lead to negligible corrections to the muon lifetime in matter. 



1. Bethe-Salpeter Analysis 



We can compute the leading bound state corrections to muonium decay using standard 
Bethe-Salpeter perturbation theory for the bound state energies ||16|| . The decay rate is 
obtained from the imaginary part of the energy (r = — 2Im£'). The leading effect comes 
from the /i + — > e + uu — > ji + contribution £ e +j, p (p) to the muon's self energy (see figs. |] and 
[1(a)). Since the decay products are typically highly relativistic, we can Taylor expand the 



self energy about the mass-shell momentum in the usual fashion |29 

p 

Im Y, e+vP (p) = — ^- + (p ■ 7 - m M ) 5Z e+uD H (Al) 

where V ^ is the free muon decay rate and m M is the muon pole mass. In muonium, the 
second term in this expansion gives the leading decay rate correction due to the bound state 



7 



reduction of the e + vi> phase space. The shift in the bound state energy due to this self 
energy is, in first order perturbation theory, 

— r — 

if> (jp e • 7 - m e ) Im Y* e + V p(pp) if> = - — ip(p e • 7 - m e)4> 

+ 5Z e+vj? V(Pe • 7 - m e) {Pn ■ 7( M ) - mjif), ( A2) 

where if) is the /i + e~ wave function, 7^) denotes Dirac matrices for the muon spinor part 
of if), p e and p M are the electron and muon momenta respectively, and integration over the 
wave function's momenta distribution is implicit. The constant 5Z e + uP is of order V ^/m^. 
The coefficient of — T M /2 in this expectation value is nearly equal to unity since the wave 



functions are normalized such that [[30 1 



V> 7(% (Pe • 7 - m e) lf) = l- (A3) 

The standard nonrelativistic expansion 



implies, therefore, that 



u{p)^ Q u{p) 2m 2 



iP(p e -l-m e )if)^(l-/3l/2) (A5) 



where /5 M is the muon velocity. This factor decreases the decay rate for muonium; it is the 
time dilation caused by the muon's motion within the atom. 

The second term in expectation value (|A2| ) is simplified by noting that 



(p e - 7 - m e ) (jv ■ 7 (M ) - m M ) V = V" (A6) 

from the Bethe-Salpeter equation, where V is the binding potential (dominated by the 
Coulomb interaction). Thus, through order a 2 , 

^(p e .7-m e )ImE e+ ^(p M )^ « (1 — /3^/2) + 5Z e + vSi (V), (A7) 

where the first correction to the free decay rate is due to time dilation, and the second 
correction is due to the contraction of final-state phase space. These two corrections are of 
relative orders a 2 m1/m 2 l and a 2 m e /m At respectively. 

A second contribution is relevant in order a 2 . This is from the imaginary part of the 
muon-photon vertex correction T p e+v9 due to the fluctuation /x + — > e + uu — > fi + - - that is, 
the correction obtained by attaching a photon to the e + in the muon self-energy diagram 
for £ e +j, p (see fig. f|(b)). Physically this corresponds to a final-state interaction between 
the outgoing positron (from the decay) and the electron (from the atom). The momentum 
transfer carried by atomic photons (~ am e ) is tiny compared to the typical e + z/z/ momenta, 
and so we can Taylor expand the vertex function to obtain the leading contribution: 

lmT p e+uP = -6Z e+uP Y + ---, (A8) 
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where the remaining terms on the right all vanish on mass shell for zero momentum transfer. 
The constant 8Z e + vV here is the same as that appearing in the self energy; this is required 
by QED gauge invariance and can be proven using Ward identities in the standard fashion. 
This vertex correction results in a perturbation given by 

-8Z e+uD ^ViP = -5Z e+vV {V) (A9) 



through order a 4 corrections (due to higher order terms in the expansion (|A8|) of the vertex 
correction). Remarkably, this final-state interaction completely cancels the part of the self- 
energy correction due to the contraction of phase space. Thus, through order a 2 ml/m 2 l the 
muonium decay rate is 

r Me = r M (1-^/2). (Aio) 

It should be emphasized that the cancelation of final-state and phase-space effects is a 
particular consequence of QED gauge invariance. Had photons been massless scalars instead 
of vector gauge bosons, for example, the final-state interaction would have vanished and the 
contraction of phase space would have been the dominant decay rate correction. 



2. NRQED and Muon Decay 

The Bethe-Salpeter formalism is awkward for analyzing muonium or other medium ef- 
fects on muon decay with greater precision. A better approach follows from the observation 
that the fluctuation /i + — > e + vv — > fi + occurs over distances, of order which are tiny 

compared with atomic length scales. Thus, the effects of such fluctuations can be modeled 
in a low-energy effective Lagrangian by local, gauge- invariant interactions. A standard non- 
relativistic effective field theory used in high-precision studies of muonium is nonrelativistic 
QED (NRQED). The corrections due to muon decay, being local, are easily included in that 
formalism as (nonunitary) corrections to the standard NRQED Lagrangian: 

£nrqed = %{iD t + - — + • • • (All) 




, ea ■ B eV • E \ 1 

+h— ir + f2— *- + ••• Wy 



O 1 J 4 Q 



Here D t and D are the temporal and spatial gauge-covariant derivatives, E and B are the 
electric and magnetic field operators, and Vv and ip e are two-component Pauli fields for 
muons and electrons. This effective theory provides a rigorous, systematic framework for 
analyzing nonrelativistic muons, including relativistic corrections as well as muonium and 
medium effects on muon decay, to whatever precision is desired. 
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The couplings ci, C2, c?i . . . in the effective Lagrangian are all dimensionless and have 
perturbative expansions in a. The q's are determined by computing muon decay, using the 
effective theory, for muons in different reference frames. They are tuned until the decay rate 
predicted by the effective theory agrees with the relativistic result 




= iyO) 1-^-^-... ; (A12) 



at tree- level c\ — — 1. The d^s incorporate the contribution from muon-electron annihila- 
tion: fie — ► vv (see Eq. (]T7D). The /j's are computed by comparing the on-shell, renormalized 
vertex function T e + uP , expanded in powers of the external momenta, with predictions from 
the effective theory. 

Two features of the effective theory deserve comment. One is that there are no terms of 
the form 

^Vje^VV' (A13) 

where A is the scalar potential in electromagnetism (that is, the Coulomb potential in an 
atom). Such a term would be required to generate a contribution such as that due to either 
phase-space contraction (second term in Eq. (|A7|) ) or final state e + e~ interactions alone. It 
cannot arise here because it would violate QED gauge invariance. This confirms the detailed 
Bethe-Salpeter analysis of the previous section which demonstrated the cancelation of such 
terms. 

The second observation is that the parameter T M appearing in the effective Lagrangian 
is precisely the free muon decay rate, to all orders in a; it is not a "bare" decay rate that 
is further renormalized by medium dependent radiative corrections. This is because the 
operator ip^ip^ in 

SC= X ^l^ (A14) 

is a conserved quantity in the nonrelativistic theory; it is the pure number operator that 
counts muons and so cannot be renormalized. Consequently 

(mm = iT* (A15) 



for any state \<p) that contains a single muon — including muonium states as well as arbi- 
trarily complicated many-electron states or other medium effects that describe muons inside 
bulk materials. This "no-renormalization theorem" implies that the only corrections to the 
free muon decay rate in such systems come from the explicit nonunitary, factorizable cor- 
rection terms in the effective Lagrangian, the terms with couplings q, di, fi ■ ■ ■ that are 
directly related to muon decay. As discussed in the text, these corrections are typically of or- 
der 1 ppb or smaller. Other effects not of this form, due to external photons or multi-electron 
excitations, cannot modify the total decay rate. 
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FIG. 1. Momentum probability distribution in muonium IS" state. 
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FIG. 2. Positron energy spectrum for free fi + decay at rest. 




FIG. 3. Self energy of the muon. Its imaginary part, computed in the external Coulombic field, 
determines the lowest order muonium width. 
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FIG. 4. Coulomb interactions in muonium decay. 
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